arXiv:1501.07576vl [math.OC] 29 Jan 2015 


Real-Time Second-Order Optimal Guidance 
Strategies for Optimizing Aircraft Performance in 
Stochastic Wind Conditions 

Kamran Turkoglu 

(.submitted and under review in Aerospace Science and Technologies) 

January 30, 2015 


Abstract 

This study presents a real-time guidance strategy for an unmanned 
aerial vehicles (UAVs) that can be used to enhance their flight endurance 
by utilizing insitu measurements of wind speeds and wind gradients. In 
these strategies, periodic adjustments are made in the airspeed and/or 
heading angle command, in level flights, for the UAV to minimize a pro¬ 
jected power requirement. In this study, UAV dynamics are described 
by a three-dimensional dynamic point-mass model. A stochastic wind 
field model has been used to analyze the effect of the wind in the pro¬ 
cess. Onboard closed-loop trajectory tracking logics that follow airspeed 
vector commands are modeled using the method of feedback lineariza¬ 
tion. To evaluate the benefits of these strategies in enhancing UAV flight 
endurance, a reference strategy is introduced in which the UAV would 
follow the optimal airspeed command in a steady level flight under zero 
wind conditions. A performance measure is defined as the average power 
consumption with respect to no wind case. Different scenarios have been 
evaluated both over a specified time interval and over different initial 
heading angles of the UAV. A relative benefit criterion is then defined as 
the percentage improvement in the performance measure of a proposed 
strategy over that of the reference strategy. Extensive numerical simula¬ 
tions are conducted to show efficiency and applicability of the proposed 
algorithms. Results demonstrate possible power savings of the proposed 
real-time guidance strategies in level flights, by utilization of wind energy. 


1 Introduction 

When nature is observed, it is possible to see that birds, are capable of taking 
advantage of wind currents not only to minimize their energy consumption, but 
also to maximize their endurance. One important aspect of this is, they do 
not hold any information and the ability to estimate the weather conditions on 
the path they fly (and/or going to fly through). All decisions are solely based 
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on with respect to local and instantaneous wind conditions. In this way, they 
optimize flight trajectory based on the local and instantaneous decisions. This 
is, a simple but, yet, inspiring mechanism to learn and incorporate into flight 
dynamics through the mechanics of flight. 

Ideally, if the regional wind information is completely known in advance 
(with the help of a pre-determined (forecasted) weather/wind maps over the 
flight region), optimal flight/trajectory planning can be used to determine flight 
paths that minimize the total power consumption over a specified time interval, 
subject to various constraints. However, the main challenge for such a pre¬ 
determined map approach is that (due to the highly complex, stochastic, coupled 
and nonlinear nature of the atmosphere) weather forecasting related prediction 
errors also propagate into the optimization routine. 

Therefore, instead of making decisions based on pre-determined weather 
maps, with this study, we propose real-time guidance strategies that will make 
local , in-situ decisions using available on-board instruments to benefit from the 
existing local wind conditions and minimize power consumption during the 
flight. The foundation of this concept had been briefly outlined in Turkoglu’s 
fundamental work mm- 

There are pioneering works in the area of UAV flights utilizing wind energies. 
The developments and flight tests of practical guidance strategies for detecting 
and utilizing thermals mm have illustrated the feasibility of these concepts. 
Patel [3] studied the effect of wind in determining optimal flight control con¬ 
ditions under the influence of atmospheric turbulence. Langelaan [4] studied 
how to exploit energy from high frequency gusts in the vertical plane for UAVs. 
In addition, Langelaan j5] presented a method for minimum energy path plan¬ 
ning in complex wind fields using a predetermined energy map. Sukkarieh et 
al. [10], TT, Hi developed a framework for an energy-based path planning that 
utilizes local wind estimations for dynamic soaring using the measurements and 
predictions from the wind patterns. Rysdyk [16] studied the problem of course 
and heading changes in significant wind conditions. McNeely m and et al. 
studied the tour planning problem for UAVs under wind conditions. McGee Q2| 
presented a study of optimal path planning using a kinematic aircraft model. 

In all of these studies, it is assumed that wind information is fully known 
over the region of flight. But in reality, wind is a stochastic process that needs to 
be addressed accordingly. One suitable approach is to devise real-time strategies 
that will benefit from the instantaneous nature of wind dynamics, on the spot, 
rather than depending on big forecasted weather maps. Thus, it is reasonable to 
adopt the idea of executing local , instantaneous maneuvers based on local/on¬ 
board measurements to utilize wind energy via the information available at that 
specific time instant t 0 . Compared with the dynamic optimization studies, this 
study presents the use and the utilization of in-situ wind measurements alone 
with no regional wind information, to optimize power consumptions. 

In this paper, optimal adjustments are made to the airspeed, heading angle 
and/or flight path angle commands to minimize a projected power consumption, 
based on the instantaneous, local wind conditions/measurements. The onboard 
feedback control system then tracks these modified (updated) commands, and 
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this process is repeated periodically throughout the entire flight. 


2 Dynamic system analysis 

When the nature of trajectory optimization problems is taken into account, 
it is possible to see that there are six main components which are vital in 
determining (and also analyzing) the flight trajectory: airspeed (V), heading 
angle (HI), flight path angle ( 7 ) and location of the aircraft, namely x, y and h. 
Once these values are known (and/or provided), trajectory planning becomes a 
relatively easy task. 

In this study, for the purpose of developing optimal guidance strategies, UAV 
flights are represented using a 3D dynamic-point-mass model, and the detailed 
structure is provided in further detail in the following sections. 

2.1 Normalized equations of motion 

In order to increase numerical efficiency and to reduce computational complex¬ 
ity, 3-D dynamic point mass equations of motion are normalized by specifying 
a characteristic air-speed V n and mass m. In this paper, characteristic nor¬ 
malization speed-V n is selected to be the maximum speed of the aircraft (i.e. 
In — Vmax)- 

Following to some algebraic manipulations [9], it is possible to obtain nor¬ 


malized equations of motion as 

v' = ^ - pV\C Do +KCl ) - siny- W / '{,(y,'F,7,5,y,ft,) ( 1 ) 

V = ^ VCl sin n- —— Wy{V,y,7,x,y,h) (2) 

cos 7 V cos 7 

— cos 'y 1— — — 

7 ' = pVC L cos/j,-+ —W^(F,1', j,x,y,h) (3) 

x' = V cos7sin + W x (x, y, h) ( 4 ) 

y' = tdcosycos^' + W y (x, y, h) ( 5 ) 

h! = ysiny-l- W h (x,y, h) (6) 

and 

Wy = W' x cos 7 sin + Wy cos 7 cos 4 / + sin 7 (7) 

W 4 = W' x cos(’F) - W' v sin(^) (8) 

W'^ = W' x sin(7) sin(l') + W' y sin(7) cos(l') — W' h cos(7) ( 9 ) 


where the functional dependencies of the wind terms are shown in parenthesis, 
for convenience. 

Since this is a constrained optimization problem, imposed constraints on 
states and controls are also expressed using normalized values, and are presented 
in the following Sections. 
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3 Problem statement 


Guidance strategies, in general, cold be grouped into three basic categories: 
action strategy, velocity strategy, and trajectory strategy. 

In this study, action strategies refer to the direct specifications of some (or 
all) of the control variables, namely zjj = (P,Cl,p) (i.e. power, lift coefficient, 
bank angle) over a certain period of time. Here, they represent open-loop control 
schemes. In comparison, velocity strategies specify some or all of the desired 
velocity components zy = ( 14 , 'F c , 7c) (be. airspeed, heading angle, flight path 
angle) over a certain time interval (At) as flight commands, where these com¬ 
mands are then followed via closed-loop tracking. Finally as its name suggests, 
trajectory strategies specify a flight trajectory of desired positions as func¬ 
tions of time over a certain time interval: x c (t),y c (t),h c (t) where t £ [to,tf]. 
These different categories of guidance strategies can be used to harvest wind 
energies of different frequencies or types, and/or be used in combination. 

A complete guidance strategy consistently needs to specify three flight com¬ 
mands at a given time, respectively along the longitudinal, lateral, and vertical 
direction. The guidance command along each direction can be from any of the 
three groups. 

Here, a high-level picture of a complete guidance strategy is provided in 
Fig-ffl 



Figure 1: High-level architecture of guidance strategies. 

In this study, it is desired to study velocity guidance strategies. The 
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combination of the listed strategies above is stil an ongoing research effort, and 
results will be reported in another studu. 


3.1 Problem formulation of real-time guidance strategies: 

Here, with the utilization of wind energy and velocity guidance strategy, it is 
desired to determine optimal adjustments in airspeed (V), heading angle ('I') 
and/or flight path angle(7) that will minimize the power consumption projected 
sometime into the future. In other words, we aim to minimize the cost func¬ 
tion expressed at the terminal cost based on local, in-situ wind measurements. 
Mathematically speaking: min zv I = P(to + At), which is subject to all appli¬ 
cable constraints. As a result, a projected power consumption expressed at the 
terminal state (i.e. at (to + At)) is used to be minimized instead of the current 
power at initial state (i.e. (to))- 

Feasible UAV flights must satisfy constraints due to UAV performance and 
operational limits. These constraints can typically be expressed as bounds on 
trajectory states and controls. For the current problem, these constraints are 


l 'i n i n ^ * . Unax 

flnin <P< -Pmax 

c Lrnln <c L < c imax 

|/^| A Mmax 

Following to the optimal solutions that are attained with respect to given 
constraints, the UAV will be directed to track optimal state commands 

V* = V 0 + AV C 


T* = T 0 + AT C 
7 * = 7o + A^ c 

For velocity guidance strategies , all corrections (AV, AT and/or A7) are the 
optimal adjustments (increments or decrements) to be determined with respect 
to local , instantaneous and on-board wind measurements. 

Once optimal adjustments in airspeed (AV), heading angle ( AT) and in 
flight path angle (A7) are obtained, it takes some finite time for the UAV to 
achieve the desired changes via closed-loop tracking. This is a phenomenon that 
has to be taken into account and it is explained in further detail in the following 
Section. 


4 Solution strategies for level flights 

4.1 Level flight strategies in presence of wind 

In level flight strategies, it is assumed that the aircraft is maintaining a level 
flight by controlling the flight path angle as zero (7 = 0). This assumption, leads 
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V(t 0 ) 


*V(t f ) 



tf - (to + At) 


Figure 2: Trajectory tracking control logic. 


to several simplifications to the presented 3D formulation, such as W' = 0. It 
also sets vertical wind component as zero (Wh = 0 ), by the nature of problem 
formulation. 

With these assumptions, at steady state conditions (i.e. V' « 0, T' ss 0 and 
7 ' w 0 ), it is possible to obtain an expression for the power at time instant to 
for level flight strategies as 

P(t 0 )n pV 3 (C Do + KC 2 L ) + VWY (10) 


In 2D level flight strategies, it is assumed that by the time (to + At), any 
commanded changes in airspeed ( V ) and heading angle (F) will be achieved via 
closed-loop tracking, as demonstrated in Figure-O 

As a result, the vehicle will be in a steady state condition: V' ~ 0 , F' « 0, 
p ft: 0 . This leads to a control command input for lift coefficient (Cl) such as 


y«0 = pVCl cos p 


cos 7 — W’ 
V 


C L 


cos 7 — W' 
pV 2 cos p 


( 11 ) 


For level flight conditions, with the assumptions of 7 ~ 0 and IF) = 0, con¬ 
trol command input for lift coefficient (Cl) in Eg. (fill) simplifies to, Cl = 
\/pV 2 cos p. Following to the derivation of control command input for lift co¬ 
efficient, using Eq. Q, power function is obtained as 


P ps VV' + pV 3 C Do + _ K , j- + VWy 
p cos^ p V 


( 12 ) 


where the projected power function (at to + A t) for level flight strategies is 
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derived as 


P(t 0 + At) « p[V 0 + AV c fC Do + j ? J AVc (13) 

+ [V 0 + AV c ]W^(t 0 + At) 

It is assumed that by the time to + At, any commanded changes in air¬ 
speed and heading angle will have been mostly achieved via closed-loop track¬ 
ing. Therefore, the vehicle will basically be at a steady state condition at the 
terminal state: V' « 0 and \x ss 0. 

Here, the W^to + At) term, which represents the main contribution on 
acceleration, needs to be expressed analytically as a function of optimal correc¬ 
tions and in the following section, the derivation of this term is described in 
further detail. 


4.2 Expression for projected wind rate in level flights 


The rate of change of any wind component can be expressed in terms of a total 
derivative as 


dW, 

'h > = Si 


dW, 


dy 


^y' + 


dW { )t . 


dh 


h! 


dW ( ) 

~dT 


where () = (x, y, h). 

In level flights, with negligible flight path angle and vertical winds (i.e. 7 = 0, 
Wh ~ 0), Wy term reduces to 


W' = 


dW x dw„ 


dy 


dx 


V sin T cos T 


0 W X - . 2 , dWy- 2t 
. V sm T H— ttt-V cos^ T 
dx dy 


+ 1 w, dw - 


dx 


W„ 


dW x dw n - 


dy 


, , T y. dW y j dW v 

+ [W x -^- + W v 


di 

dW 1 


sin T 


dx 


dy 


+ 1 cos 4- 


(14) 


Since, the only available information is the local, insitu wind measurements, 
it is assumed that the current wind gradients will stay constant m over the 
immediate neighborhood around the current position of the UAV, leading to an 
expression of Wy at (to + At), which is the projection of the wind gradient. The 
validity of the constant wind gradient assumption is strengthened even further 
with faster update rate, (i.e. smaller update time, At). 

With the constant wind gradient assumption over [to, to + At], we also have 


W x (to + At) 

Wy(t, 0 + At) 


w x . 


^-. ai+ ^ as+ ^i af 


dx 


dy 


di 


W 


Vo 


dWy A dWy A dWy , - 

V -Ax + -^AA y + -^J-A t 


dx 


dy 


di 


(15) 
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These expressions depend on (Ax, Ay), which also have a dependency on 
AV, A’F, and At through Eqs. ([4])-(J6]), and reciprocally on the wind components 
over the specified time interval. Now, in order to complete the derivation of 
the projected wind rate expression for level flight strategies, expressions for 
(Ax, Ay) will be developed in the following section. 


4.3 Expressions for position changes in level flights 


The concentration in this section will mainly be focused on developing expres¬ 
sions for (Ax, Ay). Their dependencies, explicitly, show on the increments of 
airspeed and heading angle. From Eqs. ©-©, and for 7 = 0, it is possible to 
have 


nto+At 

Ax = x(t 0 + At) — x(to) = / (V sin >F + W x ) dt 

Jio 

fto+At 


r-t o+iit 

A y = y(t 0 + At)- y(t 0 ) = (V cos 4 + W v ) 

Jt 0 


(16) 


dt 


With the assumption that both airspeed ( V) and heading angle (4) will have 
achieved their commanded values (i.e. reach steady state conditions) at the end 
of the specified time interval, and keeping in mind that wind speeds are obtained 
using Eg. (1151) (through the trapezoidal rule), the numerical integration of the 
above equations leads to 


where 


Ax — 


1 


Ay = 


Q level 

1 

Q level 


2 

At 

dWy 

dx 


, 2 \ 

Qlevel = ( ) yj 


dWy 

dy 

Bi 


2 

A t 


By 


2 

A t 


dW x 

dy 

dW x 

dx 


dW x dWo. 


dx 


dy 


dW T dW„ dW r dW„ 


dx dy dy dx 


lh 


B 2 


(17) 


(18) 


For a sufficiently small At, expression in Eg. (1151) shall always be nonzero; en¬ 
suring the existence of solutions for the position change expressions. 


4.4 Guidance algorithms for level flights 

Based on the derivations above, the power consumption at terminal state (i.e. 
at time to + At) can now be expressed as a function of the current command 
adjustments in airspeed and heading angle. Then, the problem of reducing 
future power consumptions reduces to determine AV and A*F from 

min P tf = P(i 0 + At) = I(AV, A 4; At, X 0 ) (19) 

AV.A'I' 
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subject to AVjnax and AT max incremental constraints, and the initial state 
conditions required Xq. 


5 Optimization strategies 

From the given problem formulation, it is possible to see that the nature of 
the problem is very complex, extensively coupled and highly non-linear. At this 
point, different algorithms may be used to solve the static optimization problem 
in hand. But most of the existing numerical methods heavily depend on iteration 
routines which are not desirable in real-time applications. They impose major 
drawbacks from the perspective of computation time and “convergence rate in 
allowed computation time” (which can be extremely small in some cases and 
applications). To avoid this to certain extent, in this Section, it is aimed to solve 
the static optimization problem in an analytical , and a-single-shot manner. 

For this purpose, here, yet simple but powerful gradient method will be 
presented to aid in solving the static optimization problem in hand. 

5.1 Second order optimal adjustment strategies 

With second-order gradient algorithms, the main goal is to find locally optimal 
adjustments: AG, AT and Ay using the necessary conditions for optimality. 
For this purpose, second order Taylor series approximation of projected power 
function, in the neighbourhood of initial state Xq, is taken into account 
where the analytical expressions for the local minima in the neighbourhood of in- 
situ measurements are obtained through the necessary conditions of optimality, 
General matrix representation of problem solution can be given as + 

[T 2 \ (nxn) [A]( nxl ) = 0, where for this case, corresponding components are ob¬ 
tained as 



It is possible to rewrite solution as pi]^ - \-T 2 \ (nxn) [A] (lM!l) = 0. This 
result leads to a well known least squares minimization problem of 

inf 11 y — Ax 11 
xex 

where the optimal solution is obtained as x* = (A T A) 1 A T y. 
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Therefore, in presence of favourable wind conditions, the solution to second- 
order optimal adjustment strategy (AV, AH and Ay) for the case of 2D flight 
is obtained as [A}° n p xl = [(T 2 T ) nxn (T 2 ) nxn ] _1 {~Tj) nxn (Ti)„ xl , where ( ) T and 
( ) _1 defines the transpose and inverse of a matrix, respectively. For the 2D 
level flight case, this solution reduces to calculating only AV and AH. 

This completes the derivation of optimal solution strategies and provides an 
analytical , single-shot-solution to the static optimization problem in hand, in 
presence of wind. Once optimal adjustments (AU, AH and Ay) are obtained, 
existing states are updated instantaneously 

U c * = % + AU 

V’* = H o + AH 

and optimal flight conditions are calculated. Using appropriate control strate¬ 
gies, tracking calculated optimal commands leads to a locally optimal flight 
condition that minimizes power consumption in presence of wind. Next section 
discusses such tracking control commands methodologies. 

6 Models of closed-loop tracking 

In literature, there are various (general) models of control systems to achieve 
tracking performance. One of them is the well known feedback linearization 
method. The method of feedback linearization is a widely used and a commonly 
accepted application in applied nonlinear control, and is well defined in many 
sources, as mi. In the following subsections, necessary control command inputs 
will be derived using the method of feedback linearization which will make sure 
that such desired states (U c , H 0 , y c ) will be achieved and tracked accordingly, 
to minimize power consumption in presence of wind. 

It should be noted that, with respect to the assigned mission objectives, the 
aircraft could be guided follow velocity guidance strategies (I4,H c ,y c ), trajectory 
strategies {x c ,yc,h c ) or the specific combination of both strategies. From that 
perspective, it is necessary to derive such control commands, which will make 
sure the corresponding strategy is executed as desired. Thus, in the following 
sections, necessary control commands for velocity guidance strategies will be 
derived in detail. 

6.1 Control commands for velocity guidance strategy (V a i r , 

7 ) = 

The method of feedback linearization simply defines a transformation from the 
nonlinear system in hand into an equivalent linear system through change of 
variables and a suitable control input. Here, we would like to acknowledge the 
fact that effects of plant and modelling uncertainties will have an important im¬ 
pact on the outcome of such methodology. Therefore, how those plant/modeling 
uncertainties are handled is an important matter. It is a topic which should be 
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c 


L (vel) 


\J \Wy - V cosyA'^'E - ^ c )] 2 + [cos 7 - W\ - VI< 1 ( 7 - 7c ) 


pV 2 


( 23 ) 


investigated in further detail. However, it is a topic which is out of the scope 
of the study presented in this paper and is left aside for future research. As a 
result, such uncertainty analysis is not included in this paper. 

In follows, the procedure for the derivation of control inputs ( power-P , lift 
coefficient-Ci and bank angle-/i) is presented in terms of normalized EoMs. 

If the derivation is started with first order approximation of air-speed track¬ 
ing, it is possible to define the closed-loop power law as 


= -K v (V - V c ) + pV 2 ( C Do + KC 2 Lv ) + sin7 + Wy (21) 

Similarly the closed loop control command for bank angle, (p) is derived through 
bank angle expressions as 


, _ Wj - VcosWKviV - W e ) 

an Hvel) cog 7 _ ^ _ y Kl ( 7 - 7c ) 


( 22 ) 


where the closed loop control logic for lift coefficient (Cl) is obtained as given 
in Ea. (l2dl) . 

In the given derivations above, the normalized feedback gains ( Ky , A'^, AT 7 ) 
can be tailored to reflect typical closed-loop UAV control characteristics. It is 
also good to keep in mind that, even if control commands are derived for a 
general 3D case, for 2D level flight strategies they simplify further with f »0. 


7 Simulation evaluation and results for level flight 
strategies 

7.1 Guidance algorithm parameters 

Performance of the proposed level flight guidance strategies strongly depend 
on the bounds (constraints) imposed on the adjustment parameters, such as 
AVjnax, AVmin, AlVax, A4' m i n . It is crucial to come up with ranges of their 
appropriate values for constrained optimization problems. 

Actual values used in the guidance strategies can be smaller than the given 
bounds, but in this simulational study, we select AV m ax < 5 [ft/sec], AVmin = 

—AI4iax, Af max < 30°, and AlVin = —A^max- It is also important to note 
that those values are specific to an UAV and is subjected to vary from aircraft 
to aircraft. 
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7.2 Evaluation criterion 


Because of the reason that guidance strategies provided in this study are intro¬ 
duced to save power in UAV flights, a basic performance measure is defined as 
the average power consumption over a specified time interval, 

1 -- - 

P = ■=- / TVdt (24) 

tf Jo 

where tf is the time period of evaluation (i.e. flight/simulation time). 

In the numerical results that are presented in the following sections, sampling 
time (integration step size) is selected as 50[if*], (i.e. 0.02[sec]). 

Furthermore, different initial heading angles of the UAV result in different 
relative angles with respect to the wind field, thus affecting the overall result of 
power saving benefit. In order to filter out these differences caused by different 
initial headings of UAV flights, and to have an estimate of global benefit, the 
above basic performance measure is further averaged within a set of different 
initial heading angles over [0,360°]. Here, a generic case of Ad'o = 5° increments 
are considered. The mean of basic average power consumption over different 
initial heading conditions is defined as the measure of the performance 

1 v* _ 

^ avg = iVvp ^ (25) 

2=1 

where each i corresponds to a different initial heading angle, and Ny is the 
number of different initial heading angles used. 

To evaluate the proposed guidance strategies, following two scenarios are 
considered: 

• Scenario-0: is the reference strategy that aims to follow the reference 
airspeed and a constant heading angle command set at the initial heading 
angle. This provides the reference average power consumption, P , in case 
of no wind. 

• Scenario-1: is the case where both airspeed and heading angle commands 
are adjusted periodically based on the current wind measurements. The 
resulting average power consumption is P 3 . 

7.3 Benefit criterion 

To assess the outcomes of proposed strategies, the following benefit criterion 
is introduced as a relative measure of potential fuel savings of the proposed 
guidance strategies over the reference strategy 

fr= P ° vg p 7 0 PlvS i= 1,2,3,4 (26) 

•avg 
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7.4 Simulation parameters 

In simulations, UAV parameters similar to those of the ScanEagle UAV are 
used, where characteristic values, taken from m Here, the power available is 
assumed to be able to vary instantaneously (%100 change in lsec for propeller- 
driven engines). In addition, available control rate bounds are taken into account 
as (Ptnax = 6-216, (% 100 change in lsec ),{C' Lmax = 1.865 [1/sec]), {fj,' max = 

1.085 (10 [deg/sec]). 

Furthermore, the minimum airspeed constraint is selected to be close to the 
stall speed, whereas the maximum speed is selected to be close to a typical 
cruise speed, with allowances for transient dynamics in both cases. In terms 
of normalized quantities, V min = 1.0/^pGV max , Umax = where C Lct is a 

typical cruise lift coefficient, corresponding to an angle of attack of 3° — 5°. 

7.5 Simulation results 

For the conducted simulations in this section, a level flight at altitude of 15,000[ft] 
is assumed and a total flight time of 20 [min\ = 1200 [sec] with a simulation step 
size of 0.02[sec] = 50[Pz] is taken into account. Normalization velocity, V n , is 
taken as the maximum flight speed of the UAV ( V n = 134.5 [ft/sec]). Wind mag¬ 
nitude (W m ) and wind direction ('k UJ ), obtained from the stochastic wind model 
formulation, taken from El, is used. Wind frequency, throughout these equa¬ 
tions is assumed to be constant and equal to ui w = 0.05 [rad/ft] « 9A[deg/m]. 
It is also assumed that periodic optimal adjustments are applied once in every 
4[sec], which is the same as the measurement update time (i.e. At = 4[sec]). 

7.5.1 Scenario-0: Constant airspeed, constant heading angle 

In this scenario, it is desired to fly with optimal airspeed (determined for no 
wind case) and constant heading angle, which is the initial heading angle given 
in the simulation settings. With this scenario, if the mission is specified to 
fly with optimal airspeed and initial heading angle (say 90° (with respect to 
true north)), UAV dynamics will execute these commands throughout the en¬ 
tire flight and will determine the power consumption during this specific flight 
routine. Obtained power consumption value is the reference power setting, Po, 
and will be compared with other outcomes to measure the efficiency of proposed 
strategies in given flight conditions. 

7.5.2 Scenario-1: Optimal airspeed, optimal heading angle in pres¬ 
ence of wind 

It is intuitive that combined affects of both (airspeed and heading) adjustments 
will help to harvest the maximum amount of benefit from the proposed guidance 
strategy and local wind conditions. 

For this purpose, extensive simulations have been conducted, and obtained 
results are given in Fig.®. 
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Figure 3: Relative benefit versus wind frequency response of Scenario-0 and 
Scenario-1. 


Fig-El compares the relative benefits of the two aforementioned strategies in 
which periodic adjustments are made in the commands of airspeed alone (dashed 
line), and both airspeed and heading (solid line), respectively, over the reference 
strategy. It is possible to see from Fig-[3]that the maximum benefit is obtained 
for the case where adjustments both in airspeed and heading angle are applied. 
In this scenario, overall power savings goes up to %8.5 in terms of total power 
consumption. 

It is worthwhile to note that for different wind and atmospheric conditions, 
flight path (due to heading angle adjustments) will also change. Eventually, 
while we seek to benefit from the wind currents and aim to extract energy 
from the wind, this may potentially bring the UAV to a completely undesired 
location, totally unrelated to the assigned flight mission. Thus, it is important to 
impose boundary control and restrict the flight region to complete the guidance 
strategies. 

This concept is currently under investigation and will be reported in another 
study, in further details. 


8 Conclusions 

This paper presents real-time UAV guidance strategies that utilize wind energy 
to improve flight endurance and minimize power/fucl consumption. In these 
strategies, airspeed and/or heading angle commands are periodically adjusted 
based on the insitu measurements of local wind components. It has been shown 
that using local, instantaneous wind measurements, without the knowledge of 
the wind field that the UAV is flying through, it is possible to benefit from 
the wind energy, greatly enhance the performance and increase flight endurance 
of the UAV. Throughout this research effort, UAV has been modeled using 
3D point-mass equations. Corresponding performance and practical constraints 
has been introduced to mimic a realistic flight of an UAV. A stochastic wind 
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model has been taken into account to simulate the true nature of the wind. 
UAV flights were formulated as a non-linear optimization problem and a cost 
function has been introduced to model power characteristics at terminal state 
as a terminal state cost function, which minimizes overall power consumption. 
This optimization problem has been solved as a single shot optimization, in 
real-time. Second-order gradient algorithms are used to find local, optimal 
solutions (adjustments) that will minimize the power with respect to taken local, 
instantaneous wind measurements. Extensive simulation results show that it is 
possible to obtain power savings up to %10 with respect to the flight scenario, 
with no wind. 

The proposed strategies offer improvements over the constant airspeed refer¬ 
ence strategy in terms of average power consumptions even in a constant wind 
field. These benefits initially increase as the spatial frequency of the wind field 
gradually increases, but reaches a peak at certain frequencies and then start to 
decrease beyond these frequencies. 
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